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1 Introduction 



The holographic correspondence between (super)gravity theories on AdS 
spaces and (super) conformal field theories has passed many tests and gener- 
ated much insight into the strong coupling behavior of field theory. For pure 
superconformal theories one can calculate many correlation functions from 5- 
dimensional supergravity, a procedure which is greatly facilitated by the high 
5*0(4, 2) symmetry of the bulk AdS$ geometry. A number of domain wall so- 
lutions of supergravity have also been found - 5-dimensional geometries with 
the symmetry of the 4-dimensional Poincare group- and general arguments 
show that these are dual to d — 4 superconformal theories deformed either by 
addition of relevant operators to the Lagrangian or by vacuum expectation 
values of such operators. For a Lagrangian deformation, conformal symmetry 
is explicitly broken and one expects that the trace of the stress tensor Ty and 
the perturbing operator O are related by T\{x) = (3oO{x), where (3o is the 
beta function for the operator O. For deformation by vacuum expectation 
value, conformal symmetry is spontaneously broken and one expects that 
T\ = 0. Correlation functions of the stress tensor should thus be a useful 
probe of the physics of holographic RG flows, yet there is a history of diffi- 
culty, briefly reviewed below, in attempts to calculate correlation functions 
which display the expected physics from supergravity. 

The purpose of this paper is to outline how these difficulties are resolved 
by the formalism of holographic renormalization previously developed for 
AdS bulk geometry and various linear perturbations in AdS [p]]. This for- 
malism embodies the duality between UV divergences in the boundary field 
theory and IR divergences of the on-shell supergravity action. The IR diver- 
gences are determined by near boundary analysis of the classical supergravity 
equations of motion. The bulk theory is then regularized by adding coun- 
terterms, expressed as integrals of local expressions in the fields at boundary, 
to cancel these divergences. This procedure yields a finite renormalized ac- 
tion which is a functional of boundary data for the bulk metric and other 
bulk fields. In the AdS /C FT correspondence the boundary data are sources 
for dual operators and UV finite field theory correlation functions can be 
obtained by functional differentiation with respect to these sources. These 
correlators obey field theory Ward identities including conformal anomalies. 

Near boundary analysis is sufficient to resolve all divergences, but leaves 
certain non-leading coefficients in the asymptotic expansion of solutions un- 
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determined. These coefficients contain full information on the behavior of 
correlation functions at separated points in which much of the physics re- 
sides. To find them one needs a full solution of the equations of motion, 
usually specified uniquely by requiring that the solution vanish in the deep 
interior of the bulk geometry. A full solution of the nonlinear equations with 
general boundary data is far too difficult, but one can linearize about the 
background domain wall and, if fortunate, find explicit fluctuations which 
play the role of the bulk-to-boundary propagator in AdS geometries. Two- 
point correlation functions^ can be found using these fluctuations. A simple 
cutoff method in momentum space 0, [3] is sufficient in many cases, but 
fails to give a full account of the stress tensor correlators in RG flows. We 
show how to determine these correlators from the renormalized action. Our 
formalism also simplifies the calculation of one-point functions. 

Attention was first called to the stress tensor correlators in M. The 
linearized bulk field equations couple metric fluctuations to those of the scalar 
fields which flow in the domain wall backgrounds. It was shown how to 
decouple the equations for a general domain wall, and explicit fluctuations 
were found for two flows - the GPPZ flow || describing a supersymmetric 
mass deformation of M = 4 SYM theory and a Coulomb branch (CB) flow || 
[7j dual to spontaneous breaking of the gauge and conformal symmetry of Af = 
4 SYM by a vacuum expectation value (vev) for the lowest chiral primary 
operator (CPO), bilinear in the fundamental scalars. Even with fluctuations 
known, the simple cutoff method failed to give physically reasonable two- 
point correlators. Progress was made in ||. The main difference from [|J 
was a different choice of gauge for the bulk fluctuations which enabled the 
calculation of the correlation function (T£(x)Tij(y)) for the GPPZ flow. A 
gauge invariant formulation of the problem is clearly desirable, and this is 
incorporated in the method presented below. 

There is other literature on fluctuations and correlation functions in the 
flows just described, mostly for uncoupled fluctuations such as transverse 
components of the bulk metric || [U], [TT[] , or transverse bulk vectors dual to 
conserved currents JT2| |. The fluctuations of a number of other bulk fields, 



both coupled and uncoupled, were found in [13], in which the implications of 



4 In principle rt-point correlators, for n > 3 can be calculated through Witten diagrams, 
but the integrals encountered are difficult in the reduced symmetry of domain walls. To 
our knowledge they have never been attempted, and we respect this tradition. 
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super symmetry and the multiplet structure of the fluctuations were empha- 
sized. Correlators in Coulomb branch solutions of 10-dimensional Type IIB 
super gravity were discussed in |TJ, [TJ|, [TB . 



The flows we deal with are supersymmetric. They are dual to boundary 
theories with unbroken SUSY, and (Tij) = is thus required. In field theory 
calculation with a supersymmetric regulator automatically gives (T^) = 0, 
but this relation would generically fail with a non-supersymmetric regula- 
tor. We find an analogous situation in our work in supergravity. The near 
boundary analysis used in the holographic renormalization procedure does 
not distinguish between SU SY and noia-SU SY solutions of the field equa- 
tions, so manifest supersymmetry is not guaranteed. On the other hand 
the renormalized action is ambiguous to the extent that a restricted class of 
finite local counterterms can be added. Requiring that the on-shell action 
evaluated on a bulk supersymmetric solution vanishes selects such a finite 
counterterm and (T^) = is then automatic. 

Near boundary analysis of the field equations is straightforward but quite 
complicated and differs in detail from case to case depending on the dimen- 
sion of the scalar field in the flow and its potential in the bulk action. We 
therefore try to be clear on the logical steps involved in the application of 
the holographic renormalization method to RG flows. But we simply present 
the asymptotic solutions of the field equations in an Appendix, with details 



of the procedure to be explained later [17 



The plan of the paper is as follows. In Sec 2 we review the construction 
of supersymmetric domain walls in D = 5. Our two examples can be lifted 
[§, [7J |TJj to solutions of D = 10 Type IIB supergravity. We argue that 



a Weyl transformation should be made to a frame in which the AdS^ scale 
becomes L 2 = a' (i.e. L — 1 in string units). This facilates comparison of 
energy scales in boundary and bulk and resolves a puzzle concerning the size 
of the mass gap in the Coulomb branch solution of || [7j. We also outline 
how the holographic formalism leads to the definition of the vevs (T^) and 
(O) as functions of the sources. Higher point correlation functions can be 
obtained from these quantities, and we derive the important Ward and trace 
identities which they satisfy. In Sec 3 we explain how to use holographic 
renormalization to determine (T^) and (O) explicitly in the two flows we 
study. Physical vevs are obtained at this stage of the program. In Sec 4 we 
give a general gauge invariant treatment of the linear fluctuation equations. 
In Sec 5 we present solutions of these equations and use our formalism to ob- 
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tain two-point correlation functions. Their physical properties are discussed. 



2 Holographic RG-flows 

Holographic RG-flows are described by domain-wall spacetimes with scalar 
fields turned on. For one such active scalar with canonical kinetic term, the 
relevant part of the supergravity action isQ 

S = ( d'xVG \-R + ^-GTdfidvQ + V{$)] - \l (2.1) 
JM L4 2 J 2 JdM 

where K is the trace of the second fundamental form. We work in Euclidean 
signature. 

We are interested in domain wall solutions of the equations of motion 
(EOM's) resulting from (|2.1| ) that preserve 4D Poincare invariance. For one 
convenient choice of radial coordinate, they take the form 

ds 2 = e 2A{r) 5 ij dx l dx j + dr 2 
$ = $(r) (2.2) 

The specific flows we consider are supersymmetric. From Killing spinor 
conditions in the bulk supergravity theory one can deduce the first order flow 



equations [I5| 



-£2->.,, -gWw (,3, 

where VK($) is the superpotential. The potential V($) is expressed in terms 
of W by 

V(*) = ±(d*W) 2 -±W 3 . (2.4) 

It is usually straightforward to solve the first order equations and any such 
solution automatically satisfies the second order equations of (|2.1|) for domain 
walls (but not conversely). 

Supersymmetry guarantees that the domain wall solution is stable. How- 
ever, the first order system can be derived from the requirement of gravita- 
tional stability of asymptotically AdS geometries even when the action (|2 . 1| ) 

5 Our curvature conventions are as follows R^p 17 = dpT lJ p a + T^x^T^n^ — /j, <-> v and 
Rp,v — Rp\i> X - They differ by an overall sign from the conventions in 
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is not the truncation of a bulk supergravity theory. A generalized positive 



energy argument [20], was used in [21], and it was shown that V(<&) must 



have the form ( |2.4| ) when there is a single scalar fieldQ The argument in 
[ f2l][ implies that there is a "superpotential" W such that the critical point of 
V(<&) associated with the AdS geometry is also a critical point of W.f\ In the 
AdS/CFT correspondence, positivity of energy about a given AdS critical 
point is mapped into unitarity of the corresponding CFT. It follows that in 
all cases the dual CFT is unitary, the potential can be written as in ( |2. 4| ) 
23|| . When ( |2.4|) holds a simple BPS analysis [^, |22| of the domain wall 



action yields the flow equations (|2~3|). 

We assume that has a stationary point at $ = 0. Near this point 

iy(<3>) can be approximated by 



(2.5) 



and we assume that < \i < 4. For large positive r the domain wall solution 
is well approximated by the boundary region of an AdS§ geometry with scale 
L — 1, i.e. 

A(r) par $ « exp(— jj,r) (2.6) 

If < \l < 2, then the bulk field $ is dual to an operator O of dimension 
A = 4 — fi and the domain wall describes a relevant deformation of the CFT 
Lagrangian. If 2 < \l < 4 then the dual operator has scale dimension A = fj, 
and the bulk flow describes spontaneous breaking by the vev (O).Qf] 

In the discussion above we considered domain-wall solutions that asymp- 
tote to AdS^ spacetimes with scale L — 1, rather than the official AdS/CFT 



6 For several scalars the obvious generalization of the form (2.4) implies stability, but 
the converse is not necesarily true. 
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If one relaxes this requirement the potential can always be written in the form (2.4) 



one just views (2^) as a differential equation for W |22|. In this case, however, as the orig- 
inal critical point may not be a critical point of W, the results of about gravitational 
stability do not necesarily apply. 

8 For /X = A = 2 the scalar solution dual to a Lagragian deformation is $ rexp(— 2r). 
Such a purely bosonic mass deformation breaks SUSY. Nevertheless one can obtain 
such solutions from the first order equations. The relevant "superpotential" is given by 
W(&(r)) = -3/2 - exp(— 4r)(l/8 - r/2 + r 2 ) . This "superpotential" is not analytic in <f> 
(to obtain W r ( < i>) one needs to invert $ = rexp(— 2r)). 

9 A more detailed argument (24|] shows that operators of dimension 1 < fj, < 2 can be 
described holographically. 
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scale L = (4W 2 3 s iV) 1/4 . This is a significant point that we now discuss in 
some detail. As we review in the next sub-section, the on-shell gravitational 
action evaluated on the near-horizon solution is equated to the gauge theory- 
effective action by the AdS/CFT correspondence. The latter has a large N 
expansion of the form 

oo oo 

s„ = E n2 ~ 29 E c 9 MmNY (2.7) 

On the gravitational side, specifically II B string theory in the usual Einstein 
frame, the overall constant in front of the action is 1/16ttG^ \ where 

CT } = (2.8) 

is the D = 10 Newton's constant. So the gravitational action seems at first 
sight to have a different leading behavior. As is well known the correct 
dependence on N and ^ng s = g\ u is restored because the near-horizon D3- 
brane metric depends on these quantities. To provide a more manifest match 
to the gauge theory expansion we perform a constant Weyl rescaling so 
that the solution now becomes asymptotic to AdS^ x S 5 with AdS and sphere 
radius equal to one in string units, i.e. L 2 = a'. Newton's constant is now 
equal to 

g « o, ^ 6 i^^ 6 4») 4 (2 - 9) 

(5) 

The five dimensional Newton's constant G N is now obtained by dividing 



( |2.9| ) by the volume of the unit five-sphere, vol(S 5 ) = tt 3 . It follows that 
the overall constant in front of the five dimensional action ( |2.1|) is N 2 /2n 2 . 
The difference between (|2.8[) and (|2.9|) means that energies are measured in 



different units, and the effective a' in the two frames differ by a factor of 
y/g s N. This is the origin of the factor y/g s N present in the UV/IR relation 
derived in p6[| . 

The Weyl scaling needed is the special case (p = 3) of the rescaling used to 
reach the so-called "dual-frame" [^UQ In this frame all Dp-branesQ, not just 

10 The dual- Dp-frame is defined as the frame where the metric and the 8 — p-field 
strength couple to the dilaton the same way. 

11 Fivebranes are exceptional in that the near-horizon limit gives Minki x S 3 with a 
linear dilaton. 
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the D3-brane, have a near-horizon limit AdS p +2 x S 8 ~ p . The Weyl rescaling 
and the dilaton are not constant when p ^ 3 and because of this the solutions 
are 1/2 supersymmetric even in the near-horizon limit. It has been argued in 
[ f2~5|j that this frame is holographic. It is in this frame that the on-shell grav- 
itational action has manifestly the same leading behavior at large iV as the 
gauge theory effective action, and the supergravity variables automatically 
take into account the UV/IR relation. So in order to have manifest match- 
ing between gauge theory and supergravity computations (without worrying 
about different units) we need to do the gravitational computation in the 
dual-frame. For the case at hand this simply means that we will consider the 
solution with L = 1 (we set a' = 1 from now on), and Newton's constant is 
equal to 1/N 2 . 

2.1 Correlators and Ward Identities 

To obtain correlation functions we must go back to the second order EOM's 
of ( |2.1| ) and consider solutions with arbitrary Dirichlet boundary conditions 
for the bulk fields. As in earlier papers on holographic renormalization, we 
use a different radial variable p = exp(— 2r), where the boundary is mapped 
to p = 0. The bulk scalar satisfies the (modified) Dirichlet condition 

-p^ (0) (x) (asp^O) (2.10) 
The general bulk metric ansatz is 

ds 2 = G^dx^dx u = j^ + - p g l3 {x,p)dx i dx j (2.11) 

where the boundary metric limits to 

g lj (x,p) -> g( )ij(x) (as p -> 0) (2.12) 

In the AdS/CFT correspondence, the boundary data g^ij{x) and <j)to)(x) 
are arbitrary functions of the transverse coordinates x % and are the sources 
for the stress tensor and, respectively, the operator O in the boundary 
field theory. The correspondence is then expressed by the basic formula for 
the generating functional 

(exp(-S QFT [0(o)] - / d 4 x^g^O(x)(p i0) (x))) = exp(-S SG [g(o),<P(o)}) (2.13) 
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On the left side (...) denotes the functional integral average involving the field 
theory action Sqft[9(o)] minimally coupled to g>( ). The action Ssg[9(o), 0(o)] 
is the classical action integral of ( |2.1|) evaluated on the classical solution with 
boundary data p( ) an d <j><0)- This is the leading term in the semiclassical 
computation of the supergravity partition function. Unless confusion arises 
we will henceforth use S to denote this classical on-shell action. 

Actually S is divergent due to the behavior of the solution gij(p, x), <3>(p, x) 
near the AdS boundary. One must regularize the action, e.g. by cutting off 
the radial integration at a small value p = e, and add appropriate bound- 
ary counterterms S ct to cancel the divergences. The action S must then be 
replaced by the renormalized S^en = lini^o^reg + S ct ) in ( |2.13|) and S rcn 



becomes the generating functional of connected correlation functions. This 
process is the heart and soul of holographic renormalization and is described 
in the next section. In the rest of this section we discuss some general prop- 
erties of the correlation functions obtained from the procedure. 
By definition, the variation of the effective action is equal tc^\ 

SS Ien [g {0)ij ,(j) {0) } = J c/ 4 x v /^[i(T ij )<5^ ) + (0)50(o)] (2.14) 

The expectation values (T^ ) and (O) are functions of the sources to be com- 
puted in the next section. Multi-point correlation functions can be obtained 
by further differentiation, e.g. 

KW^W) = - 1 = /ffi'ff - (2-15) 

V#(o)(y) °9(o)yy> 

A well known feature of AdS/ CFT physics is the correspondence between 
bulk gauge symmetries and global symmetries of the boundary theory. In 
the present setting the relevant gauge symmetries are bulk diffeomorphisms. 
Some of these were used to bring the bulk metric to the form ( [2. lip . The 
remainder preserve the form ( |2.11| ), and we distinguish between diffeomor- 



phisms involving only the 4 transverse coordinates: 

Sg% } = -(V ? e + VT), 50(0) = eV^o) (2.16) 

12 Notice that the definition of (O) differs by a sign from the definition used in [Q. 
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and a one-parameter subalgebra of 5D diffeos |27| whose effect on the bound- 
ary data coincides with the Weyl transformation 

Sgfa = -2agf 0) , 50 {O ) = (A - 4)<r0 (o) (2.17) 

The counterterms needed to make the on-shell action finite are manifestly in- 
variant under the 4D diffeomorphisms, but they break the diffeos that induce 
Weyl transformations on the boundary p8| |. It follows that the diffeomor- 



phism Ward identity should hold, but the trace Ward identity is expected to 
be broken. 

It is straightforward to substitute the variations ( 2.16 )- Q2.17 ) in ( [2.1 1| ) 
and obtain the Ward and trace identities^ 

V^) = -(O)V^ ( 0) (2.18) 

(T!) = (A - 4)<j> {0) (O) + A (2.19) 

where A is the conformal anomaly. As mentioned above, it originates from 
the fact that the counterterms break part of the 5D diffeomorphisms. The 
explicit form of A will be determined below. It is useful to note here that A 
is local in the sources. One of the most elegant aspects of the holographic 
renormalization formalism is the simple emergence of Ward and trace iden- 
tities including conformal anomalies. 

The source for the stress energy tensor can be decomposed as follows, 

1 

l (0)y v (* u "(0)j) -r^' 
where 

V% }ij = 0, tif£ = 0, V7if 0)i =0. (2.21) 

All covariant derivatives are that of g(py 

Using ( |2.20| ) and ( |2. 18| )-( |2TT^ ) and partial integration we can rewrite 



Sg m = Sh7 0)ij + V {i Shf 0)j) + 6i j7 6h i0) - ViVjSHp) (2.20) 



(gT4|) as 



SS IeD = Jd A x^i K --6h^(T lJ }--6hf ;(O}V 
-i^ ( o)[V^O)V^ (0 ) + (O)V 2 (O) ] 
~6h {0) [(A - 4)0 (o) (O) +A} + 5<P m (0)) (2.22) 



13 Since (Ty) includes the scalar source term in ( ^.13| ) these identities do not quite have 
the standard field theory form. For this one must use (T^qft = (Tij) + 4>(q) (C)3(o)y- 
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The form above is quite general, but we now make two assumptions appropri- 
ate to our situation. We assume that the sources describe an x-independent 
domain wall and linear fluctuations above it. This is sufficient to study one- 
and two-point functions. In this case either V 1 {0)b = or Vj0( O )B — 0, 
where the sub-index denotes a background value. The term proportional to 
5hkJ vanishes because either the source or the vev is constant and in the 
latter case the term drops out upon partial integration. One then finds 



5S ien = I d 4 x^[--5h^(T rj ) (2.23) 

1„„ ,^_o, 1 



--5# (0) <C?>V> ( o) + g^(o)[(4 - A)0 (O) (O) -A} + 6<Pw(0)] 



Since (I?) = — 85S Tea /5h(o) the expression in fl2.23|) shows that all correlation 



functions of T\ and O can be obtained from the form of (O) and A as 
functions of the sources. This is a consequence of the trace identity ( |2.19|) . 
The two-point function of T^ has the standard representation 

(TMTui-p)) = ^ kl A(p 2 ) + ir^Bip 2 ) (2.24) 

in terms of the projection operators 7Ty = 5ij — PiPj/p 2 and 

rprp ^^f0)ij 1 / \ 1 

U ijkl = - Sh T kl = ^ik^jl + KWK jk ) ~ -KijKkl ■ (2.25) 

The transverse traceless (TT) amplitude A(p 2 ) can thus be calculated by fur- 
ther variation of the TT projection of (T^-), while B{p 2 ) = (T-(jp)Tj(—p))/9. 
The Ward identity implies that (Tij(p)0(—p)) = HijC(p 2 ) with invariant am- 
plitude C{p 2 ) = (Ti{p)0{-p)) /3. Note that (T y (p)0(-p)) is the connected 
correlator. When the background vev (O) does not vanish one must correct 
for the source term in Footnote [13] to obtain this. We will obtain these cor- 
relators, together with (0(p)0(—p)) for the CB and GPPZ flows in Sec. 
5. 



3 Holographic Renormalization 



We will use the renormalization method pioneered in |28[ and developed in 
[pi], see also for related work. In this method one regulates the divergent 
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on shell action S by restricting the p integration to p > e. Asymptotic 
solutions of the field equations with arbitrary Dirichlet boundary data (= 
field theory sources) are then obtained and used to express the divergences 
as e — > in terms of the sources. Finally one adds counterterms to cancel 
these divergences. Here we will only give the results; details will be presented 
elsewhere. 

The solutions we consider are asymptotically AdS. This means that near 
the boundary one can find coordinates such that the metric #y (x, p) in ( |2.11|) 
can be expanded in a series of the form 

g(x, p) = g(o) + 9(2) P + P 2 [9{4) + h m logp + h 2 (4) log 2 p] + ... (3.1) 

The scalar field can also be expanded in a similar fashion. Since the exact 
form of the expansion depends on the mass of the bulk scalar, it will be 
presented below for the two cases we consider. 

The next step is the near boundary analysis of the EOM's. In this process 
one subtitutes the assumed expansions into the EOM's and solves them iter- 
atively. In this way many higher order terms in the expansion are determined 
as (local) functions of the sources, but not all terms are so determined. For 
the metric, gu) is the first term which is not fully determined (although its 
trace and covariant divergence are determined). It is to be expected that 
near boundary analysis does not completely fix the solution of second order 
field equations with a Dirichlet condition on an AdS boundary. Additional 
information on the behavior in the deep interior of the space-time is required. 
The product of this phase of the procedure is an asymptotic solution of the 
EOM's in which the unspecified coefficients are simply carried as such within 
the series expansions. 

It turns out that divergences of the cutoff action are fully determined 
in terms of the sources by near boundary analysis. The divergences can be 
expressed as counter terms which are boundary integrals of local invariants 
constructed from the induced metric 7^ = -g(o)ij( x ) an d the scalar field 
$(x, e). One adds these counter terms to define the finite renormalized action 
S'ren = S Teg + S ct in which the limit e — > may be taken. 

The exact form of the counterterms depends on the specific potential of 
the scalar field. For a given potential, however, the derived counterterms are 
universal, i.e. the on-shell action will be finite for any solution of the bulk 
field equations. This is a property that any holographic renormalization 
scheme should have. 
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The renormalization procedure is ambiguous to the extent that finite local 
counter terms can be added to S ren . This corresponds to scheme dependence 
in quantum field theory. In particular, S ren need not incorporate require- 
ments of supersymmetry, since it was derived using counterterms valid for 
the most general non-supersymmetric solution of the EOM's. The particular 
requirement which need not be satisfied is that E vac = in a supersymmet- 
ric vacuum. In holography this means that S TeD = when evaluated in the 
background geometry of a solution of ( |2.3|) . We will use a supersymmetric 
scheme to fix the form of the counterterms. 

It is much easier to compute the regularized on-shell action for domain 
wall solutions (|2.2|) than for those with x-dependent boundary data. For 
solutions of Q2.3p the answer can be read off from the BPS form of the 
action in |23|, || (see (13) in |23| or (14) in p2f). One obtains 



Sbk g d, r cg = / <PxyfyW[$\ (3-2) 

where 7^ = e 2A 5ij is the induced metric at the cutoff surface p = e (or 
equivalently r = — |loge). When W / ($) is expressed as a series, as in (|2.5|) , 
the low order terms will be divergent and these must agree with similar 
terms obtained from the divergences of more general solutions. In addition, 
depending on the asymptotic behaviour in p of the scalar field, there may be 
a residual finite part. If present this must be subtracted for the scheme to 
be supersymmetric, ensuring that S rcn = in the background. The specific 
working of this mechanism will be discussed with our examples below. 

It is both interesting and helpful that one can determine some divergences 
of the general S reg by examining simple subclasses of solutions of the theory, 
in our case SUSY domain wall solutions. However, one must study more gen- 
eral solutions in order to obtain conterterms necessary to cancel divergences 
in all correlation functions. For example, there are counterterms involving 
(9j0( O ) and dig^jk, which cannot be found using x- independent domain wall 
solutions. 

3.1 Coulomb branch 

Our first example is the supergravity dual of a particular state in the Coulomb 
branch of M = 4 SYM theory. It can be obtained by turning on the 5*0(4) x 
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50(2) singlet component of the scalar field dual to CPO in the 20' of 50(6). 
We henceforth denote the active scalar field by 
The superpotential is given by 

W($) = -e~^5 _ _ e 7S = $ 2 + C($ 3 ) (3.3) 

2 2 

The domain-wall solution is given by 

v = e^, e 2A = l 2 —, d ^ = 2v 2 '\l-v) (3.4) 

1 — v dr 

The boundary is at v — 1. There is a curvature singularity at v = which 
is the origin of a disc distribution of D3-branes in Type IIB supergravity. 
See H 0, [IB], [TjJ for more details of this solution and previously studied 
correlation functions. 

The change of variables that brings the domain-wall metric to the coor- 
dinate system Q2.11J ) is given by 

1 - v = l 2 p - ~IV + 0(p 3 ) (3.5) 
In these coordinates the solution <p B , A, is given by 

<Pb = -^(-Pl 2 + W + 0(p 3 )), e 2A = 1(1 - liV + 0(p 3 )) (3.6) 
V6 6 p 18 



By inspection ( |3.2j ) evaluated on the background solution ( ^.6[ ) is diver- 
gent. It follows from the given asymptotics that only the first two terms in 
the expansion of W around $ = contribute to the IR divergences, and that 
there is no finite term. All the other terms in the expansion of W vanish 
in the limit e — > 0. Thus the counterterms needed to make the background 
action finite are given by 



5 



ct,bkgd 



/ p= /*V7 (J + * 2 ) (3-7) 



We will see that these are part of the counterterms required to make finite 
the on-shell action in general. Since there is no finite term left after the 
subtraction, the renormalized action evaluated on the background solution 
equals zero, as required by supersymmetry. 
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To obtain the general form of the counterterms we note that the asymp- 
totic expansion for a bulk scalar field of AdS mass m 2 = —4 dual to an 
operator of conformal dimension A = 2 reads, 



$(x,p) =plogp(0 (o) + (2) p + p log p^ (2) 



P(<t>(0) 



(3.8) 



Inserting the asymptotic expansions in the bulk field equations one finds 
that the coefficients shown in ( |3.1| ) are uniquely determined in terms of g(o) 
and 0(o) except for gu)ij that is only partially determined and 0( O ) that is 
undetermined. In particular, only Trg( 4 ) and V l g^)ij are determined. The 
exact expressions are given in the appendix. The coefficients gu)ii and 0( O ) 
are related to the holographic one-point functions in the presence of sources 
[|I|] as we will derive shortly. 

Knowledge of the asymptotic solution allows one to evaluate the regular- 
ized action and obtain the divergences. These can be cancelled by adding 
the following covariant counterterms 



S, 



ct 



d XyFj 



---R-— loge( Rijffi - -R 2 ) 
2 8 32 6 v 3 3 ' 



+ 



$ 2 (x,e) + 



$ 2 (x, e) 
loge 



(3.9) 



where 7 is the induced metric at p = e. All curvatures are of the induced 
metric. Notice that the term <3> 2 / log e is divergent in the limit e — > with 0( O ) 
fixed (but equal to zero when 0( O ) = which is the case for the SUSY domain 
wall solutions). Thus the set of counterterms contains (|3.7f) and more. 
The renormalized action is equal to 



>Sren — lim (Sreg + S, 



ct) 



(3.10) 



where S Teg is the on-shell action in ( [2.1| ) regulated by restricting the range of 
integration to p > e. 

The expectation value of the operator dual to $ is given by 



(O) 



1 5S ri 




SS r , 



(3.11) 



By straightforward computation of the variational derivative [17 one obtains 

(O) = 20 ( o). (3.12) 
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One can similarly compute [TJ^ the expectation value of the stress-energy 
tensor. By definition 




After some computation one finds that all infinities cancel and the finite part 
is equal to 

(Tij) = g(A)ij + g [Tr gf 2) - (Tr # (2) ) 2 ] ff(0)<i - ~{g 2 2) )ij (3.14) 

1 1 ~ 2 3 

+ ^ 9(2)ij Tr 5f (2) + -(0( O) - 30 ( o)0(o) )</(o)y + ^(f> 2 (o)9(o)ij + -jh(*)v 

The last two terms can be cancelled by adding finite local counterterms to 
the action. The last term is proportional to the stress energy tensor derived 
from an action given by the gravitational conformal anomaly ( |3.22|) [0J . The 



next-to-last term is proportional to the stress energy tensor derived from an 
action equal to the matter conformal anomaly (|3.23|) . 

The computation of (O) and (Ty) was carried out in a specific coordinate 
system. One may wonder whether these results are sensitive to our specific 
choice. From the derivation of (O) and (T^) as functional derivatives of S Tcn it 
follows that they transform as tensors up to the contribution of the conformal 
anomaly. The contribution of the conformal anomaly to the transformation 
rules is also straightfoward to obtain in all generality |], Here we only 



discuss how (O) and (Ty) transform under a constant rescaling of the p 
variable. This is of particular interest because the much-discussed association 
of the bulk radial coordinate with energy scale in the boundary field theory 
suggests that the transformation 

p = p p 2 , x 1 = x 1 ' p (3.15) 

introduces the RG scale p. This transformation is an isometry of AdS space- 
time (with a metric given by ( |2.11| ) with g(o)ij = Sij and all other grmj = 0). 



Under this transformation most coefficients pick up overall factors of p ac- 
cording to their dimension, but there are also non-trivial transformations due 
to the logarithms in the asympotic solutions. One obtains 



(o) 



V) = p 2 (j)(Q)(x'p), g[ 0) (x') = g {0 )(x'p), g[ 2) (x') = p 2 g (2 )(x' p) 
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0' (o) (x') = n [0(o)(x'/x) + log// 0( O )(a//i)], (3.16) 

2 
3 ( 



2 1 

£(4)(y) = /i 4 [c/(4) + log/i 2 (/i (4) - -0(0)0(0)^(0)) - (log^ 2 ) 2 -0( O )£(o)](a:V) 



It follows that 

(O(x'))' = /i 2 ((O(xV)) + log ^20(o) (x'fi)) (3.17) 
(T ij (x , )) / = /i 4 ((Tijix'fi)) + log fi 2 [h m - <p 2 {0) g m ]{x'n)) 

It is satisfying to see that the new terms can be obtained from the following 
local finite counterterm, 

SM = J d 4 x^logfi 2 ^A (3.18) 

where A = A gra v + A sca i is the conformal anomaly given in (|3.22|) and (|3.23|) . 
This implies the transformation ( |3.15| ) only adds contact terms to correlation 
functions, and scales the momenta by 

From the expressions in ( |3.12j ) and ( |3.14| ) we find the vevs with all sources 
equal to zero are given by 

2 iV 2 

(0)b = "Ty^^ (Th)b = (3.19) 

The term N 2 /2n 2 is the overall constant in front of the action^ discussed in 
section 2. The vevs show that the solution describes a state in the moduli 
space of vacua, as promised. As one might have expected the size of the vev, 
I 2 , is set by the size of the symmetry breaking effect; the radius of the disk 
distribution of D3-branes. 

It is straightforward to use the solution of the bulk field equations given 
in appendix A.l to show that 



V l (T tJ ) = -(O)V i (o) (3.20) 

1 1 
16 v 3 3 



(Tt) = -20 {O) (O) + -( J R^-- J R 2 ) + 20 2 O) (3.21) 



14 Here and henceforth we adopt the policy of including this factor only when final results 
for correlation functions are given. 
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i.e. (Tij) correctly satisfies the diffeomorphism and trace Ward identities. 
The last two terms in (|3.21|) are what we called A in ( |2.19|) . The second 



A grav = ^(RijR ij - ^R 2 ) (3.22) 



term 



is the holographic gravitational conformal anomaly |28| and the last term 

A sca i = 2^ 0) (3.23) 

is the conformal anomaly due to matter j31] . The coefficients in both of them 
are known not to renormalize, and indeed we obtain the correct value. The 
Ward identities and the anomalies are important checks of the intermediate 
computations and of the consistency of the formalism. 

3.2 GPPZ flow 

Our second example is the supergravity dual of a M = 1 supersymmetry 
preserving mass deformation of Af = 4 SYM theory ||. We will consider 
only the simplest case in which the active scalar field is one of the two SO (3) 
singlet, dimension A = 3 scalars studied in ||. Specifically we consider 
the field called m, here renamed $, which is dual to a chiral fermion mass 
operator, and we do not treat a more general flow involving m and the second 
scalar a. 

The superpotential reads 



1 + cosh —= 



,V3, 

The domain-wall solution is given by 



|-I$ 2 -l$ 4 + 0($ 5 ) (3.24) 



V3 1 + Vl-u 2A u du 

VB = —log- 7^=, e =- , 7 = 2 1- 4 3.25 

2 1 — vl — u l — u dr 

The u variable is related to the p variable by u = 1 — p. Near the boundary 
the solution has the expansion 



¥b = P 



1/2 



[V3 + p^= + 0(p% e 2A = -(l-p). (3.26) 



^3'-^'^ p 
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By inspection one finds that the on-shell action ( p.2[ ) evaluated on this back- 
ground is divergent, and that only the first two terms in the expansion of W 
contribute to the IR divergences. To cancel them we add the counterterms 



S. 



ct,bkgd 



d xJ^j 



- + 



(3.27) 



With the addition of these terms the on-shell action is finite, but not zero 
because the $ 4 term has a finite limit. The subtraction of finite terms cor- 
responds to a choice of scheme. We proceed by subtracting the finite term 
so that the on-shell value of the action is zero when evaluated on the back- 
ground, as required by supersymmetry. In other words, we supplement the 
counterterms in ( |3.27| ) by the finite counterterm 



S, 



ct,fin 



(3.28) 



We now proceed to obtain the general counterterms. For a scalar of AdS 
mass m 2 = —3, dual to an operator of conformal dimension A = 3, the 
asymptotic expansion is 



(3.29) 



$(x,p) = P 1/2 [0(O) +P(<f>(2) +logp^ (2) ) + ...] 

The asymptotic solution can be found in the appendix. The coefficient 
is undetermined and only the trace and divergence of gu\ are determined. 
The counterterms needed in order to cancel all divergences are given by 



(3.30) 



loge 



32 



1 -r 2 ) + l -{$n^ + l -m 2 ) 



where D 7 is the Laplacian of 7. These counterterms contain ( 3.27 ), as they 
should. We further supplement them with the finite counterterm in (|3.28|) . 
This corresponds to choosing a supersymmetric scheme. 
The renormalized action is defined 



Sren — hmfSVeg + Set + ^' 



ct,fin I 



The holographic one-point functions are equal to 
1 5S rm ,. 1 1 1 SS rm \ 



(O) 



-2(0(2) +^(2)) + -0( O ) 



(3.31) 



(3.32) 
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and 

( T ij) = 9m + ^[Tr gf 2) - (Tr g {2 )) 2 ] 9(0^ ~ |(^ ( 2 2) )ii (3.33) 
1 11 

+ -g(2)ijTr 9(2) + [0(O)(0(2) - -^(2)) - ^(V0( O )) 2 ]c/(O)ij 

+^V i (O )V i (O ) + ^(^ + ^) 

where + is the stress energy tensor of the the action equal to the 
integral of the trace anomaly A = A gra v + A sca i (see appendix A. 2). The 
gravitational anomaly is given in ( 3.22|) and the matter anomaly in ( |3.3S| ) 
below. Although these terms are scheme dependent and can be omitted, we 
prefer to work with (|3.33|) keeping all the terms. 

Under the transformation ( |3.15| ) the solution transforms as follows 

0(o) W) = Mo)(^V) 

0(2) ( x> ) = ^[ < P{2){.x 'n) + log/i 2 ^ { 2)(x'/i)] (3.34) 

9(o) i x> ) = 9(0) {x'v) , 9(2) W) = V 2 9(2) {x'n) 

g[ 4) (x') = fi 4 [g (4) + log/i 2 (i^. + -7$ - 0(o)^(2)0(o))](sV) 



It follows that 



(O(x'))' = -21og^( 2 )(x / M)) (3.35) 

(T^x')) = f ( (^(xV)} + I log/i 2 (?S + TfMx'v) 



2 

Exactly as in the CB case, the new terms can be obtained from a local finite 
counterterm equal to A/2, i.e. ( |3.18| ) but with A = A gra v + A sca i (A gra v and 
Ascai given in ( p.22| ) and ( |3.39| ), respectively). 

Evaluating these expressions on the background one obtains 

{Tij) B = 0, (0) B = (3.36) 

The first of these was guaranteed because we added S ct fi n to enforce su- 
persymmetry, but this same addition was crucial to obtain (O)b = 0. The 
latter is required by the physical interpretation that the GPPZ flow at a = 
corresponds to a Lagrangian deformation of M = 4 SYM without vev. 
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Using the solution of the bulk field equations given in appendix A. 2 one 
shows that 



V i (T ij ) = -(0)V j 



(o) (O) + l(^-iit;2) 



N2 1 



[(V0 (o) ) 2 --i?^ o) ] 



(3.37) 
(3.38) 



These are the expected Ward identities. The second term in fl3.38| ), is the 
holographic gravitational conformal anomaly ( 3.22 ) [28| and the last term, 



A, 



cal 



-^[W(o)) 2 - \b4U 



(3.39) 



is a conformal anomaly due to matter [[Jl|]. The integrated anomaly should 
itself be conformal invariant, and indeed (|3.39 ) is equal to the Lagrangian 
of a conformal scalar. The coefficients are also the ones dictated by non- 
renormalization theorems. 



4 Linearized analysis around domain-walls 

As discussed in the Introduction we must go beyond the near boundary 
analysis to fix the undetermined coefficients and obtain correlation functions. 
A solution of the non-linear EOM's is beyond reach, but for our purpose 
of computing two-point functions it is sufficient to consider the linearized 
problem. 



We look for fluctuations around the solution (A(r),(pB(r)) of (|2.3|) . The 
background plus linear fluctuations is described by 

ds 2 = e 2A(r) [6 i:i + hijix, r)}dx l dx j + (1 + h rr )dr 2 
$ = (f B (r)+(f(x,r) (4.1) 

where h^, h rr and (p are considered infinitesimal. This choice does not com- 
pletely fix the bulk diffeomorpisms. One can perform the one-parameter 
family of 'gauge transformations' 

r = r' + e r (r',x'), x i = x H + e\r', x') (4.2) 

with 



/ dr'e-^d^ (4.3) 

Jr 
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where we only display the fluctuation-independent part of e\ These dif- 
feos are related to those which induce the Weyl transformation ( 2.171 ) of the 



sources. The gauge choice is also left invariant by the linearization of the 4g? 
diffeomorphisms in ( 2.16Q . 



We have seen in section || that the - part of g(o)ij drops out from the 
variation of the effective action. Equivalently the Ad diffeos can be fixed to 
set = 0, see ( |2.16|) . We thus decompose the metric fluctuation as 



hij (x, r) = hjj(x, r) + 5 i:j -h(x, r) - didjH(x, r) (4.4) 

Near the boundary the fluctuations hfj(x, r), h(x, r) and H(x,r) admit an 
expansion similar to g(x,p). The field theory sources are the leading p- 
independent parts. Analogously the scalar fluctuation <p has the expansion 
given in ( |5T5| ) or ( |3.29j ) with source ^( )(x)Q. 



The equation for the transverse traceless modes decouples from the equa- 
tions for (ip, h, H, h rr ). The coupled graviton-scalar field equations in the 
axial gauge where h rr = were derived in [|J, and we now include h rr . The 
fluctuation equations are 

[d 2 r +4A'd r + e~ 2A n]f(x,r) = 0, h% = hf 0)ij f(x,r) (4.5) 
16 

ti = -— ip'tp + AA'hrr (4.6) 
o 

H" + AA'H' - \e- 2A h - h rr e~ 2A = (4.7) 

2A'H' = \e~ 2A h + \\w v ($ - W w {p - \w v h„) (4.8) 

where is transverse, traceless, and independent of r. 

The fluctuations (if, h, H, h rr ) transform under the 'gauge transforma- 
tions' in ( [4.21 ). One can form the following gauge invariant combinations 

R = h rr -2d r (^-), h + — — 0, H'-—e- 2A q> (4.9) 
\W V ) ' 3 W/' W v * K J 



A few words about our notation: as in ( PD and (TO|), ipB and (f have a p-expansion. 



The corresponding coefficients will be denoted by fB(2k) an< l Ptek)! fc = 0, 1, .... Then the 
coefficients appearing in the near-boundary analysis are given by 4>(2k) — ( / 9 B(2fc) + f(2k) ■ 
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In terms of these variables the equations simplify, 



H' - -|-e- 2 ^ = 1 (2ii(W w - j»0 + #) (4.11) 
Equation ( |4.6| ) takes the form 

( h+ *Wv)' = -* WR (4.12) 



Differentiating ( [4.10[ ) leads to the second order differential equation 



R" + (2W W - AW)R! - (AW 2 - 2W^ W - — W 2 + -WW w +p 2 e~ 2A )R = 

y o 

(4.13) 

Equations Q4.10| )-( |4.12j ) are invariant under the transformations (|4.2|) , and 



this is important for the success of the present method. We will use these 
equations, but henceforth consider the theory in the axial gauge, i.e. we set 
h rr = 0. In 0] the axial gauge was imposed early. The resulting equations 
were equivalent to ( |4.iq )-( |4~l~2] ) but were processed differently. In particular 
because of the assumed conditions at the interior singularity, the boundary 
data for metric and field sources were not independent, and (at least in 
the GPPZ case) they were even non-locally related. This made it difficult 
to obtain 2-point functions, as one would have to disentangle the non-local 
relations of the sources from the true non-local 2-point function. 

5 Correlation functions 
5.1 Coulomb branch flow 

Let us now discuss the two-point functions for the trace-scalar sector. The 
differential equation for R becomes, 



R"(v) + -R\v) - ^ 1 , R(v) = (5.1) 
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where the prime indicates derivative with respect to v. The solution of this 
equation is 

R(v,p) = v a {\ - v)F{l + a, 2 + a, 2 + 2a; v) (5.2) 



where a=— ~(1 — yl + ^r). Equations (|4.10|) become 



2f — v) p z J 6vp z 

One can solve ( |5.3| ) for <^( ) (the coefficient of order p in the near-boundary 
expansion of <p(x,p)), 

1 I 2 8 

£(0) = Z/gfyo) + ^ + ~ 3^ (5 ' 5) 

where Q = V(l) + V>( 2 ) - V(l + a) - ^(2 + a) P|. It follows that 



N 2 ( 2l 2 1 16/ 2 \ 

(O) = ^{- 7 = 6 + ^g*(0) + ^(o)N(l + a) - 4^(1) - (5.6) 

Substituting this relation in fl2.23| ), using y/g(o ) = 1 + |^(o) + \p 2 H(q) 1 and 
going to momentum space one finds that (|2.23|) can be easily integrated to 

N 2 r A ( 1 p 2 
S Ten [h { o),H {0) ,(p {0) } = — J d p \-^-j=(p {0) h {0) - -j=ip {0) H {0) (5.7) 

2/2 g/2 \ 

+V(o) h-yg + V(o)(2^(l + a) - 2^(1) - ^)]J 

From here one immediately reads all two-point functions, 

N 2 161 2 
(0(p)0(-p)) = ^[40(1) - 4^(1 + a) + -±] (5.8) 

(2?(p)(5(-p)> = = 2(0) B (5.9) 

9/2 yy2 

pV<Ty(p)0(-p)> = --yg^ 2 = (0) ^ 2 (5 ' 10) 
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These correlation functions, as well as the one in ( |5.13| ), correctly approach 
the corresponding AdS correlators, including normalizations, in the limit of 
vanishing vev, i.e. I — > 0. Notice that the contact term in ( |5.9| ) and ( 5.1 (J|) 
are exactly the ones dictated by the Ward identities ( |3.21|) and ( |3.20|) . 



It remains to discuss the tranverse traceless sector. To determine the two- 
point function we need to obtain g^)ij as a function of ft-mw- The solution 
of the fluctuation equation (|4.5| ) is given by || 

f(v,p) =v a F(a,a,2 + 2a;v) (5.11) 

Expanding around v = 1 and converting to the p variable one obtains 

h m = h m (§^(1) + V>(3)] + ^ - £<Ka + 1)) (5.12) 

The first term on the right hand side yields only a contact term and can 
be omitted when computing correlators at non-coincident points. From the 
linearization of (|3.13|) we obtain the correlation function 



JV7-2 /„2/2 4 \ 

(TMTm(-p)) = IT?J— f ^i- - P -i>{a + 1) (5.13) 



Both ( |5.8| ) and ( |5.13| ) contain 0-mass poles. In( |5.13| ) the projection operator 



gives the pole term PiPjPkPi/p 2 ■ These poles have physical sign. Their con- 
tribution dominates the long-distance behavior of the Euclidean correlation 
functions and obeys reflection positivity. The poles are a manifestation of 
the expected dilaton, the Goldstone boson of spontaneously broken conformal 
symmetry. Pole residues are proportional to the vev (O)b- In addition to the 
massless poles, there is a branch cut along the positive real axis indicating a 
continuous spectrum above the mass scale 

I 

m w = — (5.14) 

a 

in exact agreement with gauge theory expectations: I /a' is the average mass 
of the TV-bosons H, and the vev that breaks the symmetry is proportional 
to P. The connected correlator (see Footnote I 



Al 2 N 2 2 
(TMOi-p)) = ^=2^,- = --<0>fl7ty (5.15) 
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has the dilaton pole but is otherwise purely local. Furthemore, 

(Ti(p)Tj{-p)) = (5.16) 

as can be seen from ( p.7| ) or from the linearization of ( p.21| ). These new results 
from supergravity are fully consistent with spontaneously broken conformal 
symmetry in the dual field theory! 



5.2 GPPZ flow 

The equation for R becomes 

R" + ' 



u{l — u) 



R' + 



2m - 1 p 2 



u(l-u) 4 



R 







(5.17) 



where the prime indicates differentiation with respect to u. The solution of 
this equation is 



R(u, p) = u(l - u)F{^ + \b,\~ \b, 3; u) 



(5.18) 



where b = \/l — p 2 . 

Using the definition of R in ( |4.9|) we find 



1 



>/3, 



<P(2) + ip(2) = <P(p) ~ V (^(0) + 32^(0)) J 



(5.19) 



where we used ipp) = P 2 (^ c P(o) + 32 fyo)) f rom ( |A.6| ) and J = 2-0(1) — ip{\ + 
|6) — -0(| — |fe) from p2[ . The linearization of (|3.32| ) then gives 



2 1 \/3 

(0) = -2(^ ( 2) + V(2)) + 7^1^(0) = ^(^(o) + ~^ h (o))J 



(5.20) 



Inserting ((9) in (|2.23| ) we obtain the action 

N 2 



S le n[h(0), H(o), <f(Q)] 



2vr 2 



2u2 
(0) 



(5.21) 
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where we have also included the overall factor N 2 /2ir 2 . This action gives all 
correlators of T\ and O. The last term is due to the anomaly. It follows that 
there is the operator relation, 

Tt = (3 (5.22) 

where (3o = —if(o)B = — v3- We will discuss this below. 
The correlation functions are equal to 



1 N 



2 



(0(p)0(-p)) = ~^P 2 J (5-23) 



(T!(j')0{-, ) ))=^p l ^J (5.2.1) 
2 Zn z 

3N 2 



(t!(p)t;(- p )) = ~^ p2{J+1) (5 - 25) 

The correlator (0(p)0(—p)) agrees with ||. The shift from J to J + 1 
in ( |5.25| ) is due to the linearization of the trace anomaly term in (|3.38| ). 
There is then a cancellation, and the correlator vanishes at the rate p 4 at low 
momentum. This is relevant to the issue of 0-mass poles discussed below. 

We now turn to the transverse traceless correlator. The solution to the 
fluctuation equation is |], M 

hl(u, p) oc (1 - u) 2 F(2 + |, 2 - |, 2; u)hj 0)ij (p) (5.26) 
which has the asymptotic expansion |£| 



h% = [l-pj + p 2 ^p 2 (4 + p 2 )(K - logp) + ■ ■ .}hf m (5.27) 



where K = ^(1) + ^(3) - ip(2 + f ) - ip(2 - f ). After careful lineariza- 
tion of (|3.33|) one finds the following contribution to the transverse traceless 
correlator 



( T ij) = 9(i)ij - Y^p 2h Jo)ij = ~ \m - j^p 2h fo)ij ( 5 - 28 ) 



where h( 2 ) and are the order p and p 2 terms in the expansion (|5.27|) . The 
first term in ( p.28| ) is the linearization of gity-j in ( |3.33| ); the second is a linear 
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correction to gu)ij arising because the background scale factor in ( |3.26| ) has 
an order p contribution; finally the last term comes from the linearization of 
the term 3T^/4 in ( |3.33| ). The desired correlation function is then 



AT 2 1 r? 
(T tJ (p)T kl (-p))TT = n S^[7^V + + j] ( 5 - 29 ) 

To this one must add the trace part obtained from ( |5.25| ) 

(Tij(p)T k i(-p)) trace = -^ij^H^p 2 ^ t ~ (5.30) 

We now discuss some of the physical aspects of our results. First we note 
that the correlators themselves are not completely new, [[|, ||, 10|. Rather, 
it is the coherent derivation via the holographic renormalization procedure 
which is to be emphasized. 

Physically one should expect that correlators become insensitive to the 
shape of the domain wall as p 2 —>■ oo and have the same limiting form 
p 2A ~ 4 \ogp as those in a pure AdS$ background. This property is satisfied by 
all CB and GPPZ correlators in this paper. With due care for conventions, 
this limit may be used to normalize all results. 

The (3q function which appears is a trivial constant. In fact since (f(o)B — 
\/3, it agrees with the classical valued in (|2.19|) . This is exactly what we 
expect here, since the GPPZ flow is dual to the deformation of Af = 4 SYM 
superpotential by the Af = 1 supersymmetric mass terms W = mEjTr ($?) 
with m = <P(q)b- The Af = 1 non-renormalization theorem implies that m is 
renormalized only through anomalous dimension of the operator it multiplies. 
However, in our case W is a protected operator in the undeformed theory, and 
there is no anomalous dimensionF 7 } The flo function is indeed just classical. 
Notice that our @q function is equal to the leading term in the p expansion 
of the "holographic" beta function of |36|, |37|] . 

In general our correlation functions have the discrete spectrum of poles 
noted in earlier discussions of the GPPZ flow. However there is a more deli- 
cate question of a 0-mass pole. One might have naively expected that terms 

16 See also Sec. 12.5 of f|. 

17 Contrary to naive expectation, the supersymmetric mass term does not involve the 
lowest scalar operator K, = £jTr(<^) in the Af = 4 Konishi multiplet fl34| , with non-zero 
anomalous dimension at weak coupling |34[ |35| that grows as (ggN) 1 ^ at strong coupling. 
Indeed the Konishi scalar is the simplest among the operators dual to string excitations 
pL p4| not captured by the supergravity limit. 
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in the invariant amplitudes A(p 2 ) and B(p 2 ) that vanish in the limit p 2 — > 
are just contact terms. However, because the correlators involve projection 
operators, terms that vanish at the rate p 2 yield zero-mass poles instead. The 
presence of such a zero-mass pole would conflict with the interpretation that 
the GPPZ flow is dual at long distance to M = 1 SYM which is a confining 
theory with mass gap. We find that both A(p 2 ) and B(p 2 ) vanish as p A at 
low momentum, and thus the correlators do not contain zero-mass poles. 

One may ask whether our results depend on the specific radial coordinate 
used, specifically the scaling p — > p?p associated with an RG transformation 
in Section 3. It was shown there that the rescaling introduces new terms in 
the stress-energy tensor which can be derived from the anomaly action which 
is local and thus produces only contact terms in correlation functions such 
as the complete (T^Tm). Since our method computes the TT and trace parts 
of (TijTki) separately one can see from hypergeometric series such as ( p. 2T\ ) 
that terms of order p 2 \ogp 2 appear. These potentially give 0-mass poles. 
However, one can check explicitly that the 0-mass poles cancel between the 
TT and trace parts of the correlator, leaving only the net p A log p 2 contact 
term from the anomaly. 



6 Conclusion 

In this paper we have developed a coherent approach to correlation functions 
of the stress tensor in holographic RG flows. The implementation of holo- 
graphic renormalization is somewhat tedious, but valid for all solutions of a 
given bulk action. The procedure gives relatively simple finite expressions 
for source dependent vevs (Ty) and (O) from which correlation functions can 
easily be obtained. Ward identities with correct anomalies are satisfied. The 
expected physics of flows describing both operator and Coulomb deforma- 
tions of M = 4 SYM is reflected in the results. There is more to be done. 
The correlation functions presumably satisfy Callan-Symanzik equations, see 
1 )58] for a recent discussion. The treatment should be extended to other im- 
portant field theory operators such as vector currents. We hope to discuss 
these questions and present more of the technical details of the procedure in 
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A Asymptotic solutions 

A.l Coulomb branch flow 

Scalar field: 



V>(2) 



0(2) 



0(2) 




□O0(O) + 



□o0(O) + 



g0(o)fl[0(o)]J - 7y^(o) - 2<W(°)) 

7T%(O)](0(O) + 0(0)) + 8(0(2) + V(2))) + 4^0(0) 



(A.l) 



Metric: 




30 



-4 </(2)y Tr g {2 ) ~ g (20 (o) + (o) - 30 (o) 0(o))s'(o)^ 
-20 (O )V i ( o) (A.2) 
The tensor fy 4 ) is equal to |1| 



^( 4 ) ~~ -RikjiR kl + — ViVji2 - y^ 2 ^ ~ 2l RR v 



1 _o „ 1 „o 1 , 



V 2 i? + — R 2 RmR h ) 



y 96 96 32 

= ^ (A.3) 

where is the stress energy tensor derived from the action 

S a = J d 4 x^/g^Ag r av ■ (A.4) 
Agrav is the gravitaional trace anomaly given in (|3.22 ). 



A.2 GPPZ flow 

Scalar field: 



Metric: 



^(2) = -^(□o0(o ) + g^(0)) (A.5) 



9(2)ij = 2 (^j - Q R 9(o)ijj g 

= ^(4)ii + ^%0(o) - gVi0(o)Vj0 (O ) + Y^( V ^(o)) 2 5'(o)ii 

+ ^0(O)V i V j 0( O ) + ^0( O )Do0(O)fi'(O)ii 

11 1 

= /i(4)« + + 4#(o)ijO(o)no0(o) + g^0(o)) 

^2(4)y = 0, 

1 5 

Trfi-( 4 ) = -20 ( o)0(2) - -0(o) 000(0) - ^^(o) 

16 V J 9 ; 9 
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V J # (4 )ij = V J --[Tr # (2) - (Tr # (2) ) ] # (0)ij + -(g 



^9(2)ij Tr g {2) - -h 1{A)ij - -V^V.,- 



+0(O)t. 



^(V0(O)) 2 - 0(0) (0(2) +^(2)) 



-[-2(0 (2) +^ (2) )+g0(O ) ]V 4 (O) 

where Tf, is the stress energy tensor derived from the action 



(A.6) 



(A.7) 



Ascai is the matter conformal given in ( |3.39|) . is given by the same 

formula (POD as in the CB case. 
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